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Given a family of monodromy matrices To, Ti, Tk-i corresponding to 

> ■ 

integrable anisotropic vertex models of A^-i type, \i = 0,1,..., if — 1, we 

(N ' 



O 
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O 

a: 



build up a related mixed vertex model by means of glueing the lattices on 
which they are defined, in such a way that integrability property is preserved. 
The glueing process is implemented through 1-dimensional representations 
of rectangular quantum matrix algebras A (i^ x : R n S), namely, the glueing 



matrices Algebraic Bethe ansatz is applied on a pseudovacuum space with 
a selected basis and, for each elements of this basis, it yields a set of nested 

X' 

Bethe ansatz equations matching up to the ones corresponding to an A m -\ 
quasi-periodic model, with m equal to min^z K {rank£ M }. 
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I. INTRODUCTION 



There exists a deep linking between solvable two-dimensional vertex models in statistical 
mechanics and the quantum Yang-Baxter (YB) equation, where it appears as condition for 
integrability related to some basic quantities of the model |l ] 0. This linking relies on the 
existence of an underlying symmetry, the quantum group M [|J, which comes to provide 
a nice algebraic framework to study these systems. There, solutions of the YB equation 
are representations of some quadratic relations defining the quantum group structure and, 
through a process called baxterization [0], they connect with the monodromy matrices of the 
model. From an algebraic point of view, baxterization makes an ordinary quantum group 
into a YB algebra. In this way, different representations of a YB algebra lead to different 
integrable models. 

The main aim of this work is to present a glueing process of models associated to several 
YB algebras preserving the integrability of the total system. We restrict ourself to those 
YB algebras YB n coming from baxterization of the quantum groups A (R n ) §], i.e., the 
duals of U q (su n ) |, n G N. Beside algebras YB n , the glueing process also involves their 
rectangular generalizations YB njm , defined as the spectral parameter dependent versions of 
the rectangular quantum matrix algebras A (R n : R m ) 0. Here A (R n : R n ) = A (R n ) and, 
accordingly, YB n n = YB ra . The process is based on the existence of algebra homomorphisms 
YB n m — > YB n p ® YB Pim , the cocomposition maps, that generalize the concept of coproduct 
in a bialgebra. Such maps can be used to build up representations of a given YB algebra YB n 
as a product of representations of another algebras YB m , m ^ n, in an analogous way as 
the standard coproduct is used for building up usual tensor representations. More precisely, 
given families T M and £ M , fi = 0,1, ...,K — 1, of representations of YB„ M and YB n ljnM , 
respectively, cocomposition maps ensure operator J mix = £ ® T <8> ••• <8> (k-i <%> 7~k-i is a 
representation of YB njf _ 1 (symbol (g) will be defined in next section). If each T M defines the 
monodromy matrix of a given vertex model, we say T mix is that of the mixed model with 
glueing matrices 
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We shall see this procedure is compatible with the algebraic Bethe ansatz method for 
solving these models, in the sense there exists a set of pseudovacuum vectors with respect 
to which these technics can be applied. Moreover, we show a set of nested Bethe ansatz 
equations identical to the ones corresponding to an A m _i quasi-periodic model, with m equal 
to min^ e i K {rankC^}, is related to each one of these vectors. 



This work is organized as follows: in section II, we review some well known facts on the 
connection between YB algebras and integrable vertex models; in section III, we describe the 
glueing process and the glueing matrices as 1-dimensional representations of the rectangular 
YB algebras; finally in section IV, we prove complete integrability of mixed models, showing 
that diagonalization of mixed transfer matrices reduce to solve nested Bethe equations of a 
family of A-type vertex models. 

II. YANG-BAXTER ALGEBRAS AND INTEGRABLE VERTEX MODELS 

To start with, we describe briefly the connection between two-dimensional vertex models 
and A n _\ type solutions of the YB equation. 

Let us consider the class of YB operators or constant .R-matrices 

qS k S kl , a = b; 

5 k 5 l b +(q-l/q) 5 l a 5 k , a<b; 1 < a,b, k,l < n; (1) 



5 k 5L a>b: 



J a u bi 

q G C\ {0, 1}, related to the standard Hopf algebra deformations of the simple Lie algebras 
A n _i, i.e., the quantum groups U q (su n ) and A(R n ), n G N. Baxterization process yields 
the spectral parameter dependent versions R n (x) = xR n — PR~ x Pjx of each R n , with 
P kl = 8\8 k the permutation matrix and x G C. Then, for every N G N a related integrable 
(inhomogeneous) lattice model @] is defined by a monodromy matrix T = T < - n,Ar - ) (x; a) with 
entries (sum over repeated indices convention is assumed) 

Tl = R b J(x/a )®Rll(x/a 1 )®...®R b bN _ l {x/a N - 1 ), l<a,b<n, (2) 



being a = (ckq, a^-i) a vector of C . Operators (x) : C n — > C n are entries of a matrix 



R = R fx) such that 



KM 



= [R n (x)] J in the canonical basis of C n . Compact notation 
T = R ® ... ® R, where <S> denotes matrix multiplication between consecutive factors, will 
be used. Note that T^- 1 ) = R n = R. These models are the anisotropic analogous of the 
v4 n _i invariant vertex models with periodic boundary conditions. Quasi-periodic versions 
(see again) are given by elements T in the symmetry group of R n (x), i.e., T G GL (n) 
and [R n (x) , T <g) T] = 0. Related monodromy matrices read T T = T • T. Equation (H) 



defines operators : (C r 



\<8>N 



(C 



that describe the statistical weights assigned to 



each vertex configuration in a given row of the lattice, graphically, 
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If lattice has N' rows, the partition function is Z = trace (i N 'j , being t = ^ a the transfer 
matrix. On the other hand, the operators (x; a), as it is well known, give a representation 
of the YB algebra related to R n (x). This algebra, which we shall indicate YB n , is generated 
by indeterminates T{ (x), 1 < i, j < n; x G C, subject to relations 



[Rn [x/yt] Tj; (x) Tf (y) = Tj (y) T* (x) (x/y)]^ ; 1 < i, j, r, s < n. 



(3) 



These relations entail the formal integrability of the system. In fact, by taking the trace, 
one gets [t (x) , t (y)] = for all x, y e C, i.e., the transfer matrix is a generating function of 
conserved quantities. Beside this, the model is effectively solved by means of algebraic Bethe 
|1(J and references therein), where the central ingredient is the existence 



ansatz (see 



of an eigenstate uj G (C™)®^ of each entry T^ (and consequently of the transfer matrix t), 
such that T^ uj = for all a ^ b and a > 2, and T\ uj ^ coo, Vc G C, for all b > 2. For latter 
convenience, let us express T in the block form 



T 



A B, 

Do. 



; 1 < i,j < n - 1, 
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i.e., define A = T\, Bj = T{ + , Q = Tj +1 and Dij = T^. Then, uj is an eigenstate of A and of 
each diagonal entry Da, fulfilling Q uj = and uj = for i ^ j. A vector satisfying these 



conditions is called pseudovacuum vector. On the other hand, since CUJ 
for all j, each Bj (x) plays the role of a creation operator. Applying them repeatedly on uj 
(varying j from 1 to n — 1 and x satisfying the so-called Bethe equations) we generate new 
eigenstates for the transfer matrix, namely the Bethe vectors, giving a priori a complete 
set of eigenstates for t. In such a case we say the system is exactly solvable or completely 
integrable (see [[0]] and refs. therein). Nevertheless, sometimes not only a vector but a 
pseudovacuum subspace (cf. ||13|| ) is needed in order to insure complete integrability. 
This will be our case. 

Since each YB n is a bialgebra, with coalgebra structure 

A : j{(x) » Tf(x) ® T{ (x) , e : l\{x) - S{, (4) 

for every couple of monodromy matrices T*™'^ and T( n < ps > as above we have another one, 

j(n,N + P) = j(n,N) g, T (n,P) j wlth entrles jH^N+P) = jc (n,N) ^ jb (n,P) ? 

giving again a representation of YB n . Furthermore, if uj and <fi are the pseudovacuums of 
~[{n,N) anc j j(n,p)^ then k;(g)0 defines a pseudovacuum for T( ra,Ar+p ). Consequently the enlarged 
model, or the glueing of T( n,JV ) and T^ n,p \ is also integrable. In particular, thermodynamic 
limit N — > oo preserves integrability. But, can we glue models which give representations of 
different YB algebras, e.g., YB„ and YB m with n ^ m, and such that a pseudovacuum exists 
for the resulting model? The aim of this paper is to answer last question. More precisely, 
we build up from a family {T M : \i G Zk} of pure models, i.e., T^ = T^'^, n M , 6 M, a 
mixing of them by means of glueing the lattices on which they are defined, in such a way 
that resulting mixed model can be solved by means of algebraic Bethe ansatz technics. 
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III. THE GLUEING PROCESS 



For any pair (R n ,R m ) of matrices ([I]), there exist an associated quadratic algebra 
A (Rn : R m ). They are called rectangular quantum matrix algebras 0. There are also 
parameter dependent versions, the algebras YB n m , generated by indeterminates T{ (x), 
l<i<n,l<j<m and x G C, and defined by the quadratic relations 

[Rn lx/ y )]» 11 (x) Tf (y) = T) (y) T* (s) [R m {x/y)) r ' , (5) 

1 < i> J ; < 71, 1 < r, s < m. Obviously, YB nn = YB n . In the same way as for the constant 
case [14j], there exist homomorphisms 



A p : YB n , m -> YB n p ® YB p , m ; n,m,pE N; (6) 

inherited from the cocomposition notion of the internal coHom objects, enjoying the coas- 
sociativity property (A p ® id) A r = (id ® A r ) A p |Tj|. In the n = m = p cases, these reduce 
to the usual comultiplication maps [see Eq. (Q)]. In particular, we have morphisms 

YB m > YB m n (g) YB n ® YB H)m g) YB m 

for all n, m. Now, consider pure monodromy matrices T*™'^ and T*" 1 '^ related to YB n 
and YB m , and representations A and (3 of YB m n and YB n m , respectively, where A and 
(3 denote rectangular matrices whose coefficients are representative of the corresponding 
generator algebra elements. Mentioned morphism implies A ® T^ n ' N ^ <8> (3 <E> T^ m ' p ^ gives a 
representation of YB m . As we do not want to add new degrees of freedom others than the 
related to the original models T^ n ' N ' and T^ m ' p \ we ask A and (3 to be constant (i.e., spectral 
parameter independent) 1-dimensional representations. In this case A ® T^ n,N > <g> (3 <g> T^ m,p ' 
gives an operator on (C n ) 0Ar <g> (C m ) 0P , which we shall call the glueing of T^ n,N ^ and T ( - m,p ^ 
through matrices A and (3. It is worth remarking that this is not the glueing operation 
defined in ||. Physically, A and (3 define vertices with statistical weights 

\ b and Pf, 1 < a, d < m and 1 < b, c < n. 



In general, for a family of pure monodromy matrices as described above, we can define a 
mixing of them, namely 



T mix = A <8> T ® Ai ® Ti 



(7) 



where each A M is a constant 1-dimensional representation of the rectangular YB algebra 
YB n n (mod K). Graphically, 
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being 1^ and J M multi-indices for spaces (C" - ' 1 ) M on which each acts. Since the quadratic 
relations ([5]) and the cocomposition maps (|]), one may see that ~[ mix provides a representa- 
tion of YB„ K _j. This is a direct consequence of the algebra map 



<g> YB no <g> YB no , ni ® ... ® YB nK „ 2jnjf „ 1 ® YB„ K _,. 



(8) 



Of course, these representations are highly reducible in general, as we shall see later. 



A. Constant one dimensional representations of YB n m 

Representations A^ appearing in ([/]) match exactly with 1-dimensional representations 
of A (R n : R m ), i.e., rectangular matrices A G Mat [n x m] in C such that 

[Rnt] K A? = A; \1 [R m ]Z ; 1 < i, j < n, 1 < r, s < m. (9) 

We are considering the same parameter q ^ 0, 1 for all involved i?-matrices. Otherwise, the 
only solution to ([]) is the trivial one. Using explicit form of R n given in ([I]), last equation 
is equivalent to 
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* 3 



XT \ r . = 0, 1 < r < m, 1 < i < j < n, 
0, 1 < r < s < m, 1 <i <n, 
0, 1 < r < s < m, 1 < j < i < n. 
First and second lines imply coefficients of A in a given column and row, respectively, are 
null except for almost one of them. Last line says, if ^ then all coefficients A^ with 
i < a, b < j , and with a < i, j < b, are null. Thus, each solution A of (§) is a diagonal matrix 
to which columns and rows of zeros were added. From that it follows immediately the set of 
solutions for all m, n form a semigroupoid, or a category, generated by the abelian groups 
T> n of invertible n x n diagonal matrices, and also by matrices af e Mat [n x (n + 1)] and 
<9™ G Mat [(n + 1) xn], i — 1, n + 1, n 6 N, given by 

i«(i-l)x(i-l) 0(i-l)x(n-i+2) 
x(n-j+2) 



<7„ 







(rt-i)x(j-l) -^(n-j)x(n-i+2) 



<9 n 



^<%-l)x(i-l) 0(i-l)xl 0(i-l)x(n-i) 
0(n-i+2)x(i-l) 0(n-i+2)xl Id(n-i+2) X (n-i) 



has the form 



(10) 



being 0„ )m the n x m null matrix. In fact, a general solution of Eq. 

A = d T r l dl D al ... a^ 1 e Mat [nxm]; 
a, 6 > 0, m — a = n — b = k > 0, 
with ii < ... < z a < m, ji < ... < j;, < n, and D 6 P^. If a (resp. 6) is equal to zero, then 
factors of type a (resp. d) do not appear. Such a solution has k non null entries equal to 
the diagonal elements of D, a number a of null columns in positions i a , and b null rows 
in positions ji, ...,jb- Note that rank A = k. 

Matrices cr™ and d™, which give solutions to (|9|) for m = n + 1 and n = m+1, respectively, 
are related each other by matrix transposition, i.e., d™ = (cr™) 4 , and enjoy relations 



a n ~ l a n 



a n+l Qn+1 = Qu a n_^ 

Id, 

n+1 f)n+l c\ n _n 

°i — u i-l a j ) 



a 



i < j; 
i < j; 
i =j; 
i>j + l. 
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In spite of these relations, they do not define the simplicial category. Note that, for instance, 
cr" dj +1 7^ Id. Nevertheless we name A = V n ,mgN ^n,m the category formed out by them. 

On the other hand, as it is well known, the group of diagonal matrices T> n defines precisely 
the symmetry group of R n , given by matrices D £ GL (n) such that [R n , D <g> D] = 0. 
Moreover, they are also the symmetry group of R n (x) or R n (x,y) = R n (x/y). Let us 
mention that, when i?„ (x,y) is changed by a similarity transformation Q (x) <g> Q (y) such 
that Q l k (x) = 5 l k x 2l / n , the group enlarges to T> n x Z n 0. This is why systems related to 
such i?-matrices were called Z n -symmetric vertex models |16| . 

Elements D £ D n give rise to multiparametric solutions (id ® D)~ R n (D £g> id) of the 



YB equation [17], and related twist transformations of original quantum groups |L4| [18 



Associated integrable models, which differ from the original ones by a twisting of the bound- 



ary conditions, were described in ||19|| . We shall see latter that also in mixed models the role 
of matrices D is to make a twist on the boundary conditions. 

The commutation relations between elements of T> — VneN ^Vi and A can be written 

D^ = ^Df, d?D = Dfd?, 
<- 1 D = D- oT\ D dr 1 = dr 1 Dr, 
being Df = diag (d±, df_i, 1, d^ d n ) and = diag (di, d^i, d i+1 , d n ) whenever 
D = diag (d±, d n ) £ T> n . It is worth mentioning that A can also be expressed as a product 
X = C D', where £ £ A is obtained from A by taking D = Id, and D' £ T> n is the result of 
passing D to the right, through matrices ex's, using commutation rules (|TT|) . 



B. Equivalent forms for a mixed monodromy matrix 

From results of last section, it is clear that any mixed model has a monodromy matrix 

A • R ® Ai • Ri ® ... ® Ajv-i • Rjv-i, (12) 

with = R™" for some related dimension n v , and X u = ( U D U , where ( u is in A n „_ 1>ni , and 
D u in T> nv . Here v £ Zjv- Equation (^) corresponds to the case in which there exist K 
numbers iV^, \i £ Z^, giving a partition of and such that 
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R M , = - = Ra/.+^-i = R" M - and \ u = ld Uy for v ± M , M K ^, (13) 

being M = and M M = J2a=o N a for 1 < /i < K — 1. Furthermore, Eq. (|12|) can be brought 
to an equivalent form 

T m » = Co • Ro ® Ci ■ Ri ® ® Civ-i • Riv-i T, (14) 
where T is an element of V nN1 . To see that, let us define matrices D^' £ P njfc , £ Zjv, by 



n(fc) 



Jd nfc , < fe < z/, 



A/, = z/, 

and for each > v by the solution of D v ( v+ % ( u+2 ... Cfc — Ch-i Cf+2 ■■• Cfc ^ ■ We mean by 
'the solution' of last equation the inversible diagonal matrix that arises when passing, 
in the first member, the matrix D u to the right using (|TT|) . Then flljD and ( 14 ) are similar 
through Yi u&N ® k&N £>(*), and T = U U&N D[ N ^\ 

In other words, every mixed model is physically equivalent to a twisted version of another 
one whose corresponding matrices A„ are in the category A. The role of matrices P's is to 
implement a twisting of the boundary conditions. Accordingly, we can describe each mixing 
jmix j n t erms Q f a f am ily of elements £ A nv _ uTlv , the glueing matrices, and a diagonal 
matrix T = diag ...,r njv ^ of T> nN1 , the boundary matrix. 

Another useful expression for the monodromy matrices of these models can be given from 
the following observation. Any matrix A £ A n m of rank k (note that k < m,n) may be 
written A = PAP', where A = <9™ _1 ... <9f +1 crf +1 ... cr™ _1 £ Mat [n x m], that is, 



A 



x(m— fc) I 

j 

i 0(n—k)xk 0( n _k)x(m— fc) y 

and P, P' are appropriate permutations. More precisely, if 

A = d^ 1 ... d k h a\ ... a™' 1 , i 1 <...<i a <m,j 1 <...<j b <n 
[see Eq. fllOD l, then we can choose, for instance, P £ Mat [n] and P' £ Mat [m] to be 
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(15) 



p _ n. n, n. nri j p> _ n. n, n, 

respectively, being C rjS , r < s, the matrix that acting on the right (resp. left) makes a 
cyclic permutation sending s-th column (resp. row) to r-th one, and acts as an identity for 
the rest of columns (resp. rows). Hence, for a given family of glueing matrices we have 
( u = P v ( U P U . Introducing last expression for ( u into Equation (|T1D, and making a similarity 
transformation ® v <z% N P' v , an equivalent system 



T m« = p oCo . R() .g o (g^-Ri-gi® ... 

... ® (n-2 ■ • Qn-2 <2) Cjv-i • R/V-l -FW-1 ^ > 



(16) 



where Q v = P u P u+ i and = {P' v ® P' v ) K v {P' v ® P' v ) \ follows. That is, 

f mx = {® v&N p' v ) T mx (® u&N p' v y l ■ (17) 

It can be seen for an arbitrary permutation that 

R a a (x) = R a a (x) , and R b a (x) = x 2 ^ R b a (x) for a ^ b, (18) 

where coefficients e a b takes values —1,0,1 depending on the considered permutation. We 
shall show in the next section that mixed models with 

Po = P'n-i = Idn N - x , and Q„ = Id nu ,\/v eZ N - X , (19) 

are solvable by means of algebraic Bethe ansatz technics [actually, ( |T9"D can be slightly relaxed 
and ask Q u = ld nv ^v e Z^, instead]. Furthermore, we shall see complete integrability 
implies transfer matrix obtained from ([IB]) is similar to the trace of 



jmi, = Co • Ro ® Ci • Ri ® - ® Cn-i ■ RiV-1 T. (20) 

Thus, we can solve all vertex models with glueing matrices ( u = P v ( u P' v satisfying (|TJ]) by 
solving those with J mix given in Equation (p0|). In addition, all ( u can be supposed to have 
the same rank. 
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IV. INTEGRABILITY OF MIXED VERTEX MODELS 



In order to show exact solvability of these models (or unless of a subclass of them), since 
needed commutation rules follow from map @, we must prove there exists a suitable set of 
pseudovacuum vectors for ~[ mix from which all its eigenstates and corresponding eigenvalues 
can be constructed. In other terms, using block form 



T 



("mix r\mix 



1 < i, j < riN-i - 1, 



we look for elements <3> G 'H mix = ® U £i N ( & n " which are eigenvectors of A mix and of each 
diagonal entry D™ x , such that C™ x $ = 0, and D™ x $ = for i ^ j. In this way, we build 
up recursively all eigenvalues and eigenstates by applying repeatedly operators to the 
mentioned vectors. Completeness problem will be studied separately. Of course, the smaller 
the rank of involved glueing matrices, the smaller the set of monomials in B™ x and the bigger 
the number of pseudovacuum vectors we need to construct the complete set of eigenstates. 
In the singular case for which rank^ = for some z/, we have ~T mix = and accordingly 
every vector of T-C mix is trivially a pseudovacuum vector, and no creation operator is needed 
in order to diagonalize the transfer matrix t mtx . Note in this case, operators BJ 1 ^ are null. 
Thus we can have pseudovacuum vectors which are annihilated by operators B™" and still 
be able to build up an eigenstate basis for t mix . 

We actually show exact solvability for a particular class of mixed models. Concretely, 



we concentrate ourself in monodromy matrices whose related glueings ( u satisfy Eq. (19). 



A. The pseudovacuum subspace 

Let us first consider the mixed models with monodromy matrices ~T mix defined by Equa- 
tion (pO"D, that is, each ( u = ( u is of the form (|T5|). They are a particular case of those 



with glueing matrices satisfying (0). At the end of this section the general case will be 
addressed. In order to simplify our calculations, we shall suppose 
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m = rank (o < rank ( u for all v G 'Ln- (21) 

This can be reached by a similarity transformation cyclically permuting tensor factors of the 
linear space FL mtx . Of course, m < rank^ < rain {n u ^i, n u } (mod N). Also, we suppose 
m > 0, since for m = diagonalization of T mix is immediate. Note that fl2~T|) implies 

and in particular 

o-miz v- >riJV - 1 TO mix -ra mix l\mix i \r->m— 1 pvmia; ( C ?^V\ 

t - La=l 1 a - La=l 1 a ~ A + 2^i=l D ii ■ ^ 

For a < ra we have (non sum over 6) 

where sum over each c w is in the interval 1 < c v < rank£„. 

Let us indicate by ei,...,e n the elements of the canonical basis of C n . Then Fi mix is 
spanned by vectors of the form e/ Cg> ... <8> e f N -n which can be identified with an obvious 
subset F of functions / : — > N : z/ i— »■ / y . In particular, given / G F, we denote IF the 
corresponding vector of FC mix . We shall show there exists a set of pseudovacuum vectors, on 
which algebraic Bethe ansatz will be applied, labeled by the subset Fo of functions 

/GF / Image/ C {1} U {n G N : n > ra} . (25) 

More precisely, there exist vectors $^ G Fi mtx 1 f G Fo, expanding a space 

Ho = span G : / G F } C H mi \ (26) 

namely the pseudovacuum subspace, and fulfilling 

A m« $/ = n d ni/g/ _ 1(1) G ( x / av ) $/, 

D mix $ / = r . +i d (j < m ) ? D^J = C™" $ / = 0, 

being d = l\u& N 1/G (x/a„) and G (x) = (xq — 1/qx) / (x — 1/x). In particular for the 
equal rank case, i.e., if rank£„ = m for all u, then $^ = Qf . Note that FCq C ker Q™ %x Vi. 
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(27) 



We also show 

B;" r «I> r / r<I> r if le Image/; otherwise, Bf ix ¥ = 0; (28) 

i.e., each B™* creates new states when f^ 1 (1) ^ 0. Using that we construct a set of Bethe 
vectors from each , with j from 1 to m and f^ 1 (1) 7^ 0, and generate in this way all 
eigenstates of the transfer matrix. 

To find the vectors we need some previous results. 



1. The action f~[ mix n vectors $V 

Let us evaluate the entries J b a mix on each vector Qf . From Eq. (]22|) it follows that 
jb^mix Qf _ q £ or a jj a > m. So we only consider a < m. As usual 0, we normalize 
operators = R n " b a (x/a u ) : C n " — > C n " in such a way that on the canonical basis of C n " 

^ 5 b a /G(x/a u ) e k , k ^ a, 

+ (l - ^) c S g(fe-a) {x/a v )J e b , k = a, 
being c± (x) = (q — 1/q) x ±x j [xq — 1/qx). Then, from (|24]) and the first part of fl29|), it 
follows that 

Tt m<1 0/ = r b 5 b a Uuez N l/G(x/a u ) £l f = r a 5 b dnf (30) 
if a ^ Image/, i.e., if f iv) ^ a for all ^ G Zjv- In particular 

T>™*Q/ = i/ a ^ Image /U {6}. (31) 



Also, if / G Fo and 1 < a < m, since in this case a Image/ [see (p5|) 1, we have that 
D% ix = T+\ mix Qf = r i+ i d Qf, for 1 < i < m, 

D™* fi' = TlS mx & = 0, /or 1 < i, j < m, 1 + j, (32) 

Qmix Qf = jl mix Qf = Q ; f Qr 1 < i < m> 

putting i = a — 1 in fl3"U|). Otherwise, let a a be the first integer such that / (cr a ) = a, that is, 
f [y] ^ a for all ^ < o~ a and / (a a ) = a. Let us write 
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Qf = Q9 a ®e a ®fi /a , with g a , f a : Z^, Zjv- CTa -i — »• N. (33) 

If cr a = iV — 1, we take iV" equal to 1. Then, using (p^ ) and (p5| ) again (note that a does 
not belong to Image g a ) we have 

J b a mix & = El=i kCCTa C a i Q 9a ® ei <g> ™ x 

(34) 

C<M = ^a.i IL <( 7 a V G (a/a*) ! Atf = ^ + (1 - <%) c sg(i-a) (z/ « CT J • 

Here, operators ~[ b . mix are given by the last iV — cr a — 1 factors of T mM! . From Eq. fl3"ip 
applied to ~[ b . mix Qf a y since C7 0) i 7^ for all 1 < % < rank£<j a , the non zero terms of (34) are 
those with i inside I a = {i G Image f a U {b} : i < rank^ CTa }. If b is the unique element of I a 
and b ^ Image f a , then 

T b a mix Qf = r b C a , b (Uu^ a l/G{x/a v )) W a ®e b ®nf a 

(35) 

= r b C a>b d G {x/a aa ) VL 9 * ®e b ® tt fa . 
Otherwise, suppose there exists c\ G I a fl Image / a , and let a aci be the first integer such that 
f a (cr aCl - a a ) = ci. Then f fi/ a = E, e / aci C aciji ^ ® e, ® f \ ™* ^ with 

O^i = D CUI U, a <u<a aci 1/G (ar/a„) , = ^° C1 (g) e Cl ® fi/ aci , 

iaci Image / aci U {6} : i < rankC CT(lCl } • 

A recursive process easily follows and the generic term reads 

with C aci ... Cfeii = D Cjfc)i Ua aci ... Ck _ 1 <u<a aci ... Ck l/G(x/a v ). Of course, ci G J a n Image/*, 

W-.c^ = {»6 Image r cic - c - 1 U {6} : % < rankC^...^} (36) 

and Cj G I cic 2 ...cj_i H Image / aciC2 "' Cj ' _1 for 2 < j < k. The process ends when 6 G I ci...c fc ) 
6 ^ Image / aci "' Cfc , and we choose Ck+i = b. In this case 

j^nr^- 1 = r b c aci ... Ck , b iu^..., n^'"* ® e5® n^-"*. 

In particular, writing ~r& mia; ^/ = J2 g er^ab^ 9 we h ave the given sequence of numbers 
C\, Ck G Image / defines a function g such that t{ b 7^ 0, being 
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Q9 = Q9 a $ eci $ ^9 ac i g c2 ^ ^9 ac l c 2 g ... 



c 2 ...c fc 



(37) 



Note that can be written 



Qf = fif a g e a g> fis aci g e ci ® fi 3 " 



ic 1 c 2 ...c fc _ 1 



e Cfe _ 1 <g> fi* 



tci eo.-.ci. 



e Cjk <g>fi' 



Furthermore, defining Jy s = {a a , a aci , a aci ... Cfc } and c = a, and recalling c^+i 
have [compare with Eq. ([351) 1 



(38) 
6, we 

(39) 



If # [/ 1 (a)] = + 1 and a = b, the sequence of numbers Cj = a, i — 1, k, corresponds to 
the vector Q 9 = fF. Also, Jff = f^ 1 (a) and accordingly, since D aa = 1, 

tU = r a dU v ef-Ha)G{x/a v ). (40) 
Comparing fl5?| ) and (|38]) , we see that functions g such that 7^ necessarily satisfy 



Image g U {a} = Image / U {b} 



In addition, for each element u G Image /, u 7^ a, b, function g must hold 



(41) 



(42) 



and 



# [(T 1 ^)] = # [r 1 (a)] - (1 - o ab ) , 

# fer 1 (&)] = # LT 1 (6)] + (i - * ab ) ■ 

Now, defining the classes of functions C C F, in such a way that /, g G C z/f 

Image / = Image # and # (u)] = # [g^ 1 (u) 
for all u contained in their respective images, we can write 



(43) 



(44) 
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where Cj^ is the class given by functions g satisfying (|4~1~D, ( ^2|) and (|3|). Let us note that 
C = C+a iff a = b. Then, denoting 

H x = span e H mix :/6l) for each ICF, (45) 

spaces Tic are T^; '""-invariant for all a. On the other hand, when a ^ b (since C 7^ Cj^), 
vector Qf can not be written as a linear combination of vectors fi 9 's appearing in ([0]) (they 
form a linearly independent set of vectors). That is, J b a mix Qf is not proportional to 
Also, if a ^ Image /, then C+„ = and consequently T^ mix = 0, such as follows from Eq. 
( PD|) for a ^ b. Last observations translate for operators B™ x into equations 

B m« x n / _^ cfi/; -j y-i ^ _^ 0. gmix — q otherwise. (46) 

Let us briefly study the reducibility of the action on Fi m%x of the algebra generated by 
operators T™ x . It follows from Eq. (|36| ) that, if M = max v {rank £ v }, numbers c\,...,Ck 
and Ck+\ = b must be smaller than or equal to M. This implies T b a mtx = for b > M, and 
we can restrict ourself to the a,b < M case. Also, comparing ( |37|) and (]38f) , if f (u) > M 
then g(y) — f (y). As a consequence, beside (fLT|), (fi^) and (|4*3D , condition 

g(u) = f (v) \fv e Z N such that g(v),f (y) > M (47) 

is necessary in order to have t{£ 7^ 0. Thus, defining the classes E C F as those whose 
functions satisfy pTj), it is clear that spaces are invariant under the action of T mix . It 
actually can be found smaller invariant spaces inside FLe, depending locally on the ranks of 
glueing matrices, but we will not discuss it here. 

For the equal rank case we have m = M, and accordingly the classes E are in bijection 
with elements of Fo- Thus, we can decompose Ft m%x into T mM -invariant subspaces H-Etf) 
labeled by elements of F o- In addition, by a simple inspection of coefficients fl39|), it can be 
shown the actions on Ti.E{f) and 'He{ 9 ) are equivalent provided f' 1 (1) = g^ 1 (1). Moreover, 
in the homogeneous case, namely, a v = 1 for all u, above equivalence still holds when 

# [T l (i)] = # br 1 (i)]- 
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In the following subsection we diagonalize (when possible) the operator k mix restricted to 
each He, and show its eigenvectors, when C C Fo, are precisely the pseudovacuum vectors 
we are looking for. 

2. Diagonalization of A mix and vectors 

Let us consider a class of functions C. Using Eqs. (pEDD and 
defining a fg = t{{ for / ^ g, we have that 

A m » & = a f Qf + Egec&f *fg ^, a f = r x d ELe/-^) G (x/a u ) . (48) 

Now, we are going to show there exists a total order relation between the functions of C, 
such that w.r.t. this order we can write 

A""' 0/ a f ill' ■ Y ;I . f a fg Q 9 . (49) 

In other words, operator A mix restricted Tie is represented by a triangular matrix w.r.t. the 
resulting ordered basis [recall Eq. ((4|) for X = C}. 

To see that, let us consider a function / 6 C. Assign to f u the number 1 if f u = 1 or 
if f v 7^ 1. Denote & the binary expression related to the sequence /o, /jv-i- From Eqs. 
( P7|) and (p8|) for a = 6 = 1, we see that ft* > b 9 (as real numbers) for f ^ g, since unless 
one f v — 1 were moved to the right. This implies aj g = if b* < b 9 , that is, 

A mix Qf = a f Q f + E 6g<6 / a /fl (50) 

So let us define an order < between the elements of C by saying g < f if b 9 < V , and when 
b 9 = ¥ we choose an arbitrary order. Using that and equation above, Eq. ( |49| ) follows 
immediately. 

Since eigenvalues of IK mtx are given by the numbers a/, in order to insure its diagonal- 
izability we can ask the considered model to be completely inhomogeneous, i.e., a v ^ 
for all Ijn- Then a/ ^ a g provided f^ 1 (1) ^ g^ 1 (1). Thus, eigenvalues are distinct, 

unless those related to f-th and g-th rows for which / _1 (1) = g~ l (1). But f^ 1 (1) ^ g^ 1 (1) 
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M) for a = b = 1, and 



if b f > b 9 . Therefore [see Eq. (|50])], A mix does not mix vectors related to rows with the 
same diagonal entries, and accordingly A mtx is diagonalizable. Actually, we just can insure 
fimix _ fimix ^ j g diagonalizable f or almost all values x of the spectral parameter. Note 
that for some isolated points x Q G C, we can have a/ (x ) = a g (x Q ), in spite of condition 
f^W^g-^l) holds. 

Using usual recursion formulae for diagonalizing triangular matrices, we can define for 
each subspace Tic the basis & , f G C, given by 

, \n min °, f = mtn c , , 

& = < (51) 

+ E g <fXfg® 9 , f>min c , 

with 

a g+g / (a g+ - ag), f = g + , 

Xfg = \ (52) 
[ [a fg -E g< h<fafhXhg)/ {a f - a g ) , f > g + ■ 

Here mine is the minimal / G C w.r.t. the defined order, and g + is the first element in C 
bigger than g. Equation ( p^) must be understood as a recursive formula on / for each g. 
Because [A mtx (x) , A m%x [x')\ = for all x,x' G C [that follows from commutation relations 
given in @], operators A mtx (x) can be diagonalized simultaneously. Thus numbers Xfg an d 
vectors $^ do not depend on the spectral parameter. 

Let us note diagonal entries D™ x can be diagonalized as above. But this is not enough to 
diagonalize \J nix i since operators A mix and D™ z do not commute among themselves. Never- 
theless, last operators restricted to FCo do commute, and accordingly can be simultaneously 
diagonalized. This follows from the facts that FCq C ker C™ iX and that Eq. (|3]) implies 

D mfa {x) ; fornix = {%) qrf. (y) c _ + g™z (y) Cf* (x) C + (x/y) . 

Now, let us see that vectors $^ for / G F o, given by QSTp and (0), satisfy Equations 
( p7|) and Since they are eigenvectors of A 771 ^ with eigenvalues a/, the first part of (|27|) 
follows immediately. For the second part, note is a linear combination of vectors Q 9 with 
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g inside C . Also note, if / G Fo, then the class defined by / is inside F o too. Hence, using 
Eq. ( j32|) we arrive at the wanted result. The same happens for ( |28| ) using Eq. (|46|). 

For the equal rank case, it can be shown that aj g = for all / G Fo- In fact, sets I\ Cl ...cj 
defined by fl36|) (putting a = b = 1) has 1 as the unique element, and consequently the only 
possible sequence is q = 1 for i = 1, # (1)] — 1. Such sequence correspond to the 
vector Qf . Then, the latter is an eigenvector of A mtx (without any inhomogeneity condition). 
In other words, K m%x restricted to TLq is represented by a diagonal matrix for the basis Q* , 
f G Fo, and accordingly $^ = ft* . 

To end this subsection let us say last results, valid for monodromy matrices T mix of the 
form (|20|), also holds for those given by Eq. ([16]) and satisfying condition (|19D. In fact, on 
the canonical basis ex,..., e nv of C™", using Eq. fll8|) and fl2"9|), we have that 

8 h a /G(x/a v ) e k , k ^ a, 

(5 b a + (l - $a) (x/a l/ ) 26ab c sg(b _ a ) (x/acvf) e b , k = a. 

Then, applying J b a mix to a vector $V we arrive at Eqs. (|30|) or fl34|) , depending on Image/, 
where the second term of coefficients Cj [see Eq. (|34])1 must be just changed by a factor 
(x/a a ) 2£ai . Therefore, all above results follows. In particular, all we have said for fk mix is also 
true for A m * x , and the former is diagonalizable iff so is the latter. There is a minor change 
in coefficients df g , and consequently in the linear combinations ( pH| ) that define eigenvectors 
of A mtx . Denoting the latter by $^ , and recalling Eq. (0), we conclude 

Theorem 1. Given a mixed vertex model J mix = £ • R ® ^ • R x <g> ... g> (jv-i • Rjv-i T, 
with glueing matrices = PvC,vP' v satisfying Eq. (|T9|) and (|2l|) , and assuming A mia: is 
diagonalizable (e.g., T mix is completely inhomogeneous) , it follows that vectors 

& = (®vez N Plr l $' f , / e To, 

are pseudovacuum states for T mix satisfying Eqs. ([27]) and (p8[). When rank<^, = m W, 
A™ x is diagonalizable and $ / = (g^z^) -1 fiA □ 
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All that can be rephrased in terms of our original mixed monodromy matrices, i.e., in 
the form (|7|). We just must regard them as particular cases of ( |T2| ) subject to (p~3|) . 



B. Nested Bethe equations 

Let J mix be a monodromy matrix as that given in theorem above. Thanks to the algebra 
embeddings YB n _i <^-> YB„, n > 1, which are a direct consequence of equations 

[ R n-l\ab = [ R nTab > for 1 < d,b, k, I < U - 1, 

it follows that T b a mix for a,b < m satisfy relations corresponding to the YB algebra YB m . 
Then, following for each $^ e 7io analogous technics to the ones developed in refs. J7| | 
that is, proposing as eigenstates for t" 1 ^ [see fl2Bp ] the Bethe vectors 



= B}',"" (.r,: a) ... Bf* (x r ; a) <h f . ./,. ...../V < 



m, 



and separating in the so-called wanted and unwanted terms, we arrive at a set of nested 
Bethe ansatz equations which in its recursive form are given by 

lLe/-i(i) ^ /^J 

(53) 

and 

Ayji—l (XJ 7"m) 

(*) = P ~ r \ ,J A, +1 (x) + r i+1 n^ 1 G (x^/x , {Km -2) 

where I = l,...,m — I, k = l,...,r;, and < r\ < r;_i < # (!)]• Thus, Bethe 
equations related to a given $^ e 7i , are the ones corresponding to an A m _i type 
quasi-periodic vertex model with nf = # (1)] sites per row and inhomogeneity vec- 
tor aif = (a UQ , a ui , a Un _ X J, such that z/j 6 (1) and z/j < for all z G Z n/ . When 
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# [/ 1 (1)] = we have no Bethe equations. Note in this case B™ x $^ = Vj [see Eq. ([28])]. 
For each solution 

x={xW = (4°,...,a;«):/ = l,...,m-l} 
of Equations Q53|) and (|54|), 

A/ (x; x) = d nlLi G Ax (x) + n d Uuef-HD G (x/a„) UlU G (x^/x) (55) 

gives an eigenvalue of t mix . Note that (x; x) = A 9 (x; x) if f^ 1 (1) = g^ 1 (1). This is the 
main source of degeneracy for the transfer matrix. It can be seen each A^ (x; x) differs by 
a factor EL^f-ifi) l/G(x/a v ) from the corresponding eigenvalue related to the mentioned 
A m _i model. Eigenvectors iff* (x), i.e., the Bethe vectors, can also be given recursively, but 

£rn-i\ w ith coordinates (yff^ n "' 3r i (w.r.t. the canonical basis 

of C m ~ l ) such that 

¥ (x) = - Bj" (xP; a) ... B£" a) (56) 

for 1 < I < m - 2 

and \l/ m _i = 1. Here ji,...,j n+1 < m. We are denoting by u>i the pseudovacuum for the 
pure monodromy matrix T^ m ~ l ' Tl '. Let us mention, in the Z-th level of nesting process the 
involved monodromy matrix actually is the twisting 

T (m-i,n). Th bemg T, = dia0(7i,...,7w_,), 

which also has uj\ as pseudovacuum vector. 

Summing up, we have constructed a set of eigenvectors for t mix by applying creation 
operators BJ^'s over all / G Fo- I n the following section we address the combinatorial 
completeness of that set of states. 

By last, let us say that equations (|53|) and Q53D do not depend neither on permutations 
P V -,P' V defining the glueing matrices of ~T mtx (recall conditions of theorem above), nor on 
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the set of ranks of the latter. They only depends on the minimum m = min v {rank^,} of 
that set, on the boundary matrix T, and on the inhomogeneity vector a. Hence, assuming 
complete integrability, the spectrum of the related transfer matrix t mtx , which would be 
given by the numbers A* (x;x) defined in (|55|), only depends on m, T and a. Accordingly, 

Theorem 2. Assuming complete integrability, every mixed model with glueing matrices 
satisfying Equation dl9|) is physically equivalent to one with monodromy matrix of the form 
( POD and satisfying the equal rank condition: rank<^, = m for all v. □ 



C. Combinatorial completeness 



In this section we are going to show that Eq. (|56|) (varying indices j from 1 to m, 
functions / in Fo, and x along solutions of (|53|) and ([54])) defines unless dimH mix = I\ u( zz N n v 
different vectors. That is to say, we have a set of Bethe vectors from which, a priori, a basis 
of eigenstates for the related transfer matrix can be extracted. To see that, we shall assume 
combinatorial completeness of Bethe ansatz equations related to the A n _i vertex models, 
i.e., for a model with N sites in a row we suppose there is unless (n — l) r different 
solutions for the Bethe equations corresponding to r creation operators. This has been 
shown for n = 2 (see for instance pi]), but we do not know about any similar result for 
bigger n. In our case, we would be saying for each vector 3v with / e Fo, there exists unless a 
number (m — l) r r^J of different solutions of ( |53D and ([54]) corresponding to n = r creation 



operators. Recall that n/ = # [/ 1 (1)]. Let us first see why this assumption is useful for 
our purposes. 

It is enough to analyze the case of monodromy matrices given by (p0|). The other 
cases, i.e., those given by (|HD and satisfying (|T9|), follow analogously. So let us come back 
to §IV.A.l and consider the action of operators BJ" X ' with j < m, on vectors $V with 
f E. Fo- Suppose first that rank^ = m for all v. For a = 1 and b = j + 1, sequences 
Cj = 1, % = 1, k, with 1 < k < rif define terms proportional to vectors Q 9 = fF w ', with 
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fi G f^ 1 (1), (u) = f (is) for all v \i and / M j (\i) = j + 1. That is, we change a vector 
ei by a vector e^+i in position \i G (1). They are the only possible sequences. Thus, the 
action of each B™", j = 1, ...,m — 1, on a vector Qf gives rise to a linear combination of 
rif linearly independent vectors. Existence of n/ different solutions to Eqs. (|53| ) and (pip 
for ri = 1 and for each j, is a necessary condition to obtain rif l.i. eigenstates from the 
set of Bethe vectors. Then, varying j from 1 to m — 1, we shall have, a priori, (m — 1) n/ 
l.i. eigenstates. Applying B™^ and B™ x we have (m — l) 2 vectors, each one of them having 
rif (rif — 1) /2 l.i. terms. In general, if we apply r creation operators to fl* , we have (m — 1) T 
vectors with related r^J terms. Now it becomes clear why our assumption is needed. The 
same argument can be given for the general rank case. There, when an operator BJ 1 " acts on 
we have as above the terms proportional to iV' tJ , [l G / 1 (1), together with additional 
terms given by vectors VL ha ' j with h belonging to the same class of /. Thus, the latter appears 
as terms when B™ x ' is applied to Q h . Accordingly, in order to avoid overcounting, we do not 
have to take them into account. 

Let us come back to our original problem. If combinatorial completeness holds there 
exists unless a number J2rLo ( m — l) r ("/) = (( m — -0 + l)" 7 = mW/ °^ Bethe vectors for 
each function / G Fo- Thus, since < rif < N for every / G F, the total number of Bethe 
vectors is S/eFo 771 ™ 7 = ^k=o mk Pk, being p k the number of functions / G Fo such that 
rtf = k. Let us calculate p k . It is clear that the number of functions / in Fo with the same 
pre-image / _1 (1) is 

„ fo, uEf^CL), 

II (n v -mf\ e v = (57) 
v&n 1, otherwise. 

In terms of numbers eo, ...,£jv-i) the condition rif — ^ [f^ 1 (1)] = k can be characterized 
by equality e + ... + £at_i = N — k. Then, in order to obtain p k we must sum over all 
configurations of eo, ...,en-i (s v equal to or 1), such that last condition holds, i.e., 

Pk = E eo ,..., £JV _ 1 Uuez N (n„ - raf v 5 £Q+ .„ +EN _ 1;N - k . (58) 
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Accordingly, 



But 



— Ee ,-,ejv-i IL/eZjv { n P — m T V EfeLo m?C ^e +...+ejv-i,JV-fc 
= m N T,e ,...,e N _ 1 ILeZjv - l) 



= rw* (fe - 1)° + - 1)) = u, &N *t = mr N Tl 



and consequently S/gfo mn/ = rLeZjv n v as we wan t e d to see. 



CONCLUSIONS 

From last equation we see that, under conditions of Theor. 1 and assuming complete 
integrability, T-C mtx can be decomposed into a direct sum of m n - f -dimensional spaces 7if, each 
one of them generated by the Bethe vectors related with some / inside F o- Note this sum, 
in general, is not orthogonal w.r.t. the usual scalar product in ® u ez N C' n " y . Thinking of the 
quantum spin ring related to our vertex model, whose Hamiltonian H is constructed from the 
logarithmic derivative (if there exists) of the transfer matrix, states oiTif can be interpreted 
as those of an anisotropic A m _i type spin chain with rif sites, which are localized on the 
subring Z n/ ^ f^ 1 (1) C Zjy In other words, we have decomposed a mixed spin model 
as a direct sum of A m _\ type ones with different numbers of sites and generically different 
inhomogeneities. Multiplicity of these models is given by flSTD [recall eigenvalues ( [5"5] ) only 
depend on / through f^ 1 (1)]. In connection with Theor. 2 let us say that for the equal 
rank case, since we have Tif = Ti-E{f) ( see at the end of §IV.A.l), described decomposition 
(which results orthogonal) and mentioned multiplicity are direct consequences of the facts 
that last spaces are T m ^-invariant, and that corresponding actions on spaces TCf and 7i g 
are equivalent when f^ 1 (1) = g^ 1 (1). 
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In the homogeneous case we have in addition actions onTij and 7i g are equivalent still 
when rif — n g . In other terms, for the homogeneous equal rank case we can write 'H mix as 
the orthogonal direct sum 

H mix = 0^ o CP*®W A [see (|58|) for numbers pk] , being fk some 

function with nj k = k. 

Concluding, we have presented a procedure for glueing different integrable vertex models 
in such a way that the integrability of the whole system is preserved. This procedure relies on 
some generalization of the coalgebra structure to the case of rectangular quantum matrices 
and their representations, enhancing the deep linking between these algebraic structures and 
integrability. 
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